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Q ■ Abstract 

Properties of multiparticle final states produced in the central rapidity 

K^ . region of small-x deep inelastic scattering (DIS) are discussed. It is pointed 

CN ■ out that these properties contain important information on the nature of the 

cn ' 

CN i Pomeron - used for the description of cr\*p and diffractive production in DIS. 

t^^ . It is shown that models based on universality of the Pomeron predict charged 

0\ 

■^ ' particle multiplicities in small-x DIS which are in good agreement with recent 

^ ■ 
Ph! HERA results. 
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1 Introduction 

Small-a; DIS provides a good testing ground for theoretical approaches to high- 
energy interactions of photons and hadrons. A strong increase of the structure 
function of the proton F2{x,Q'^) as x decreases observed at HERA |I|, ^ has been 
often considered as an indication of the existence of the "hard" Pomeron, predicted 
by QCD-perturbation theory [^ ^. On the other hand models based on the stan- 
dard perturbative QCD-evolution provide a good description of experimental data 
l^-P. However, these two approaches do not explain the striking difference between 
the energy dependence of total cross sections for hadronic (or 7p) interactions and 
7*p-interactions observed experimentally. In ref. a unified approach to hp, 7p 
interactions and DIS has been proposed, where the Pomeron pole is assumed to be 
universal, i.e. the same in soft and hard processes, with an intercept ap(0) ~ 1.2. 
The non-universality of the energy behaviour for hp{'jp) and 7*p-interactions is at- 
tributed to the different sizes of unitarity corrections (multipomeron exchanges) in 
these two classes of processes. It was argued that in DIS the unitarization correc- 
tions due to multipomeron exchange are much smaller than in hadronic interac- 
tions. In the latter the unitarity corrections reduce the effective Pomeron intercept 
from ap{0) ~ 1.2 to roughly 1.1. On the contrary, in DIS the "bare" intercept 
ap{0) = 1.2 is only slightly modified by the unitary corrections and provides a "sin- 
gular" (F2(x, Q^) ~ x~°'^) initial condition for perturbative QCD evolution. The 
change between an effective intercept of about 1.1 and the bare one of 1.2 occurs 
very quickly when Q^ increases from zero (at Q^ >j 1 GeV^ the effective intercept is 
already quite close to the bare one). 

A study of structure functions of the proton can hardly discriminate between 
different theoretical approaches. Investigation of diffractive production of hadrons 
by real and virtual photons provides extra tests of different theoretical models. 

Here we would like to point out that the properties of multiparticle production 
in DIS are related to the nature of the Pomeron and can be used to test theoretical 
models. The same situation takes place in high-energy hadronic interactions. A 
study of only 0")^° (s) and da}^ {s,t)/dt does not allow to discriminate between 
models with large unitarization effects |T^, |T3] and models where these effects are 



small [|1J]. However a fast increase with energy of charged particle densities in 



the central rapidity region, broad multiplicity distributions and long-range rapidity 
correlations observed in hadronic interactions clearly demonstrate the importance 
of effects related to the multipomeron exchanges \T2\ |T^ . 



We shall demonstrate below that the model based on universality of the Pomeron, 
discussed above, gives parameter free predictions for charged particle densities and 
multiplicities in the central rapidity region of DIS. These predictions are in a good 



agreement with HERA experimental data ||16 



2 Description of the model 

The CKMT-model for small-x DIS corresponds to the diagrams for elastic 
7*p-amplitude, shown in Fig. 1, where the upper blobs denote all diagrams which 
describe the coupling of the Pomeron to virtual photons. They contain in particular 

diagrams corresponding to QCD-evolution. Asymptotically (i.e. for Y = in— > 

oo), at fixed Q^, the Pomeron amplitude increases as exp(Ay) (where A = ap(0) — 1) 
with increasing energy, while the blobs in Fig. 1 have finite widths in rapidity. The 
s-channel cutting of the diagrams of Fig. 1 corresponds to the natural physical 
picture of 7*P inelastic interaction [|r^, shown schematically in Fig. 2. It illustrates 
the fact that a highly virtual photon produces a cascade of partons with virtualities 
decreasing along the chain. Systems of partons with small virtuality (of large size) 
interact strongly with the proton either inelastically (Fig. 2a) or diffractively (Fig. 
2b). This picture is similar to the models of Bjorken |TB[ and BiichmuUer [r^|. The 



inelastic production amplitude in Fig. 2a can be expressed as a sum of inelastic 
cuttings of a single Pomeron diagram (Fig. la), two Pomeron exchange diagram 
(Fig. lb), etc. Some of these cuttings are shown in Figures 3. 

In the models based on the 1/A^-expansion, the Pomeron corresponds to diagrams 
for elastic scattering with the topology of a cylinder. The cutting of Fig. 3a leads 
to a 2-chains configuration as shown in Fig. 4a and that of Fig. 3c to the 4-chains 
configuration of Fig. 4b. 

The weights with which different configurations contribute to the multiparticle 
production are given by AGK-cutting rules pO] and thus can be fixed if the relative 



contributions of rescatterings in Fig. 1 are known. 



It was already mentioned above that contributions of rescatterings in DIS are 
much smaller than in hadronic interactions and in first approximation, at availa- 
ble energies, it is possible to consider only the PP-contribution. In this case, the 
latter contribution can be determined from HERA data on large rapidity-gap events 
pT| , p2| . The ratio of the second rescattering (PP-exchange) to P-exchange is equal 
(with minus sign) to F,p /F2 (where F^ is the diffractive contribution to F2), which 
is equal to 10 ^ 15 % in the region 10"^ <x < 10"^ and 2 GeV^ < Q^ < 10^ GeV. 
For hadronic and 7p interactions, this ratio is substantially larger. 

It is known that particular multiparticle configurations are more sensitive to con- 
tributions of rescatterings than the total cross section. For example the shadowing 
contribution to the cutting of a single Pomeron shown in Fig. 3b is four times larger 
than the contribution of the diagram of Fig. lb [^. For a study of multiparticle 



production in DIS we will assume a definite model for the series of rescatterings 
shown in Fig. 1. An adequate description of hadronic interactions was obtained in 
the "quasi-eikonal" model ||23| . It is the simplest generalization of the eikonal model 



which allows one to calculate all rescatterings in terms of the double rescattering 
term. So we shall apply the same model to DIS. 



In this model the total cross section a^*p{W , Q ) can be written in the form ||24] 



a^ifiw^ Q^) = E ^?p(w^^ Q') = <yp{w\ Q')f (I ) (1) 



where k is the number of cut Pomerons {k = corresponds to diffractive produc- 
tion). The contribution of the single Pomeron exchange to cr^.p is denoted by 
ap{W\Q^) = ^(g2)exp(A0 with A = «p(0) - 1 and ^ = in^ = ln\. Here we 
have neglected the real part of the Pomeron amplitude - which has a small effect on 
multiparticle production. The function f{Z), defined as 

00 / y\n—l 

takes into account all rescatterings in terms of only one function Z = ^2!^/ K exp(A^). 
The function Z can be determined from the data on the diffractive cross section 



c^;., = ap(W^^g2) 



/ (f ) - fiz) 



(3) 



or, more precisely, from the ratio 



^D — '~^l*P 
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(4) 



The quantity {R^ + a'p^) in Z is related to the t-dependence of the diffractive 
production amplitude with R"^ = 2.2 GeV-^ and a'p = 0.25 GeV'^ [g 

The cross sections with fc-cut Pomerons (corresponding to final configurations 
with 2fc-chains) have the following form 



ct',*,{W',Q') 



kZ 



fc-l r^ 

l-exp(-Z)^ 



j=0 



2! 



i> 1 



(5) 



It follows from eqs. (0)-(^) that for small values of the ratio R^ one has Z ^ 
8 ■ R^. Experiments on diffraction dissociation at HERA show that the ratio R^ 
is practically Q^-independent and we shall consider C{Q'^) as a constant in future 
calculation. From the HERA data on the ratio R^ |^, ^, ^ we obtain C ^ 1.5 
GeV-2. 

In order to calculate multiparticle production in the chains of Fig. 4 we shall use 
the DPM-model |T^ or QGSM-model |jl3[ which give a good description of rapidity 
and multiplicity distributions in hadronic interactions. Thus we do not introduce 
any new free parameter in the calculation of multiparticle production in DIS. 

In both models the hadronic spectra of each chain are obtained from a con- 
volution of momentum distribution and fragmentation functions. The former are 



obtained from Regge intercepts [|T2| [0. In the QGSM, used in the following calcu- 
lations, the fragmentation functions into charged hadrons are given in ref. |25| with 
a normalization constant a^ = 1.14. The multiplicity distributions are computed 
assuming a Poisson distribution in clusters for the individual chains - with fixed 
rapidity positions of the chain ends. The formalism is described in detail in ref. 
fl^ . The only parameter is the average charged multiplicity of the cluster decay, 
for which we use the same value K = 1.4 as in previous works. The weights for the 
configurations, with k-cnt Pomerons (2A;-chains) are given by eq. (^. (A Monte 
Carlo code for multiparticle production in 7p collisions, based on DPM, has been 
introduced by Engel and Ranft [^). 

It is important to note that our calculations are reliable only in the central 



rapidity region and the proton fragmentation region, while the fragmentation region 
of the virtual photon can be dominated by perturbative production of jets of the 
type shown in Fig. 2. This region occupies an average rapidity ~ £n-^, out of the 
total rapidity interval ~ fra-^— . Unfortunately experiments at HERA cover mostly 
the photon fragmentation region and only a part of a central rapidity region. In the 
following comparison with experiment we will consider only the data in the central 
region of rapidity. 

3 Comparison with experiment 



Let us consider first the density of charged hadrons in the central rapidity 

dy ~ o-(*°«)(VK2,Q2) dy 

corresponding quantity for nondiffr active processes "^^ = -^ ^^^ , where a^^ 



region ^„ = „(tot)(w2 o2\ "^ h,, ■ ^^ i^ convenient also to consider the 

dy a'^o dy 

oo 

(J (to*) _ g-0 _ ^ cr«. Y)ue to AGK-cancellations [EO] the inclusive cross section in 

k=l 

the central rapidity region ^ ^ '~f^ ^^ determined only by the contribution of 
the Pomeron pole and thus behaves as (^j . For hadronic interactions, where 
rescattering effects for o"*^*"*^ or a^^ are very important, the increase of a*-*"*-' with 
energy is much slower and, as a result, ^ (as well as —f^) have a fast increase 
with energy. It was already emphasized in ref. p^ that in DIS cr*^*"*-' and a^^ 
have approximately the same ( M -behaviour as ^ (due to the small effect of the 
rescatterings) and ^ should have weaker energy dependence than for hadronic 
interactions. Likewise, multiplicity distributions in DIS are expected to be narrower 
than in pp or 7p and to show a much smaller KNO-scaling violation. 

Predictions of the model for the mean charged multiplicity of nondiffr active 
events in the central pseudorapidity domains of DIS as functions of energy iW) 
are compared with Hl-data [|I6| in Fig. 5a. Densities of particles in chains are 
taken from the model [l^. The DPM model leads to practically identical results. 
As explained above, the predictions of the model are most reliable for the smallest 
pseudorapidities (1 < r/* < 2), where they are in good agreement with experiment 
both in absolute magnitude and in the rate of increase with energy. In Fig. 5a we 
also show the results for the region 1 < 77* < 3 where the agreement between theory 
and experiment, although still reasonable, is less good. It becomes slightly worse 
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when the size of the rapidity interval increases. Note that the agreement in the range 
1 < ?7* < 3 is better at higher energies, where the distance between central rapidity 
region and fragmentation region of the virtual photon is larger. Model predictions for 
the rapidity distributions in DIS (for Q^ = 15 GeV^) at different energies are shown 
in Fig. 5b. A substantial increase of the central plateau is predicted. At energies 
W< 200 GeV this increase is mainly connected with the increase of charged particle 
densities in the Pomeron chains. The faster increase at energies y/W > 10^ GeV 
is due to an increase of the average number of chains (cutted Pomerons). At these 
extremely high energies {x < 10^^) effects of unitarization for structure functions 
are very important. Note, that as anticipated above, the energy rise of the central 
plateau in DIS is substantially smaller than the one measured in pp collisions. 

The experimental data for charged hadrons produced in nondiffr active DIS in 
different pseudorapidity intervals of the central region [|I^ are compared with our 
predictions in Figs. 6. The agreement with experiment is good in the interval 
1 < ?7* < 2 and becomes worse for larger intervals). 

Comparison of theoretical predictions for different moments of multiplicity dis- 
tributions with experimental data |TB| is given in Table I. It is important to notice 
that the values of these moments are substantially smaller than the ones measured 
in pp scattering at the same energy (see Ref. |2^). This is a clear confirmation of the 
main feature of our approach, namely that the size of the rescattering contributions 
is smaller in DIS than in pp. 

A further confirmation of this result could be obtained from the study of the 
energy dependence of the moments of the multiplicity distribution in DIS. It has 



been predicted in [|^ that in the central rapidity interval \y*\ < 1.5, the moments Cj 
in pp should have a strong energy dependence. This dependence has been confirmed 
experimentally |^ (see also |12[). Such an energy dependence is due to the increase 
with energy of the fluctuations in the number of chains. Since the multi-chain (i.e. 
multi-cut Pomeron) contributions are smaller in DIS, we expect a smaller increase of 
the moments with energy. Note that models such as JETSET and MEPS based on 
perturbative QCD predict a magnitude and energy dependence of the DIS moments 
Ci in the central rapidity 1 < r^* < 2, which is substantially different from our 
predictions |jT6|. 



Experimental studies of p_L-dependences of final liadrons at HERA |^, |31| are 
also in agreement with our assumption that in the central rapidity region [y* < 3, 
for W = 200 GeV or xp < 0.07) properties of final states are similar to the ones in 
7P or hp. Average transverse momenta in this region are much smaller than in the 
fragmentation region of the virtual photon |^ . 

Finally, let us now comment on multiparticle production in the virtual photon 
fragmentation region, where QCD perturbation theory can be applied. Though our 
model is strictly speaking not reliable in this region it still reproduces some features 
of multiplicity distributions up to pseudorapidity r^* ~ 4 though the agreement with 
experiment is worse than in the central rapidity region t]* < 2 -h 3. The model 
predicts very weak dependence of ^ on Q'^ for fixed value of W and a rather strong 
dependence on W for fixed Q^. These consequences of the model are in an agreement 
with experimental results |]T6[, which show that mean charged particle multiplicity 
is practically independent of Q^ in the interval 10 GeV^ < Q^ < 400 GeV^ for fixed 
W. This observation is not easy to reconcile with predictions based on perturbative 
QCD, - radiation of extra gluons in the final state should lead to a substantial 
increase of multiplicity as Q^ increases |3^. It may be an indication that the non- 



perturbative dynamics described above is important even in the fragmentation region 
of the virtual photon. 

4 Conclusions 

Our analysis indicates that production of hadrons in the central rapidity region 
in DIS can be understood in an approach based on the universality of the Pomeron. 
The model presented here does not contain any adjustable parameter. The size 
of the rescattering contribution has been determined from diffractive production 
data in DIS and all other parameters are fixed [Q, |13[ from high-energy hadronic 



interaction data. The model gives a quantitative description of the densities of 
charged particles in the central rapidity region and their multiplicity distributions 
at HERA energies. Possible tests of our approach in forthcoming experiments have 
also been presented. 

Our results provide a clear confirmation of the idea, introduced in previous works 



[0|, that unitarity corrections, which play a capital role in pp interactions, are 
substantially smaller in DIS. Our model can provide a quantitative way of relating 
multiparticle production in hard and soft processes. 
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Figure Captions 

Fig. 1 : 

The Pomeron-exchange diagrams for elastic 7*p-scattering amplitude. 

Fig. 2 : 

Inelastic cuttings of the diagrams of Fig. 1. a) Totally inelastic final state, 
b) Diffractive dissociation of a virtual photon. 

Fig. 3 : 

a) s-channel cutting of the single Pomeron exchange diagram of Fig. la), 
b), c) Totally inelastic cuttings of the PP-exchange diagram of Fig. lb). 

Fig. 4 : 

a) Two chains configuration corresponding to the diagram of Fig. 3a). 

b) Production of 4-chains corresponding to the diagram of Fig. 3c). 

Fig. 5 : 

a) Comparison of theoretical predictions for energy dependence of mean charged 
multiplicity of nondiffractive events in central pseudorapidity domains of DIS with 
experimental data [|16]. Lower points and curve correspond to 1 < rj* < 2 and the 
upper ones to 1 < r/* < 3. 

b) Model predictions for the rapidity distributions of charged particles in the 
central rapidity region at various energies. 

Fig. 6 : 

Multiplicity distributions of charged hadrons produced in nondiffractive events 
in central pseudorapidity domains of DIS |T^ compared to theoretical predictions, 
a) ly = 97 GeV, 1 < r]* < 2 ; h) W = 202 GeV, 1 < r]* < 2 ; c) W = 97 GeV, 
1 < r]* < 3 ; d) ly = 202 GeV, 1 < r^* < 3. 

Table Caption 

Table 1 : 

Model predictions for various moments and cumulants of the multiplicity distri- 
butions of charged particles in DIS in two central pseudorapidity intervals and at 
two different energies are compared with HI data W^. 
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<W >-- 


= 201.9 






l<r]* <2 


1 <r]* <2 


1 < r^* < 3 


1 <r]* <3 


C2 


1.69 


1.69 ± 0.02 ± 0.03 


1.48 


1.44 ± 0.01 ± 0.02 


c^ 


3.81 


3.75 ± 0.12 ± 0.20 


2.83 


2.62 ± 0.09 ± 0.09 


C4 


10.73 


10.02 ± 0.65 ± 1.22 
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5.71 ± 0.47 ± 0.34 


R2 


1.34 


1.32 ± 0.02 ± 0.03 


1.29 


1.26 ± 0.01 ± 0.03 


R3 


2.26 


2.11 ± 0.09 ± 0.13 


2.06 


1.88 ± 0.08 ± 0.11 


Ks 
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0.147 ± 0.043 ± 0.101 
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0.108 ± 0.041 ± 0.036 
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1 < r^* < 2 
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1.71 ± 0.02 ± 0.02 


1.48 
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3.85 


3.87 ±0.16 ±0.18 
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4.83 ± 0.19 ± 0.64 
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1.19 ± 0.01 ± 0.03 


R3 


2.11 


2.11 ± 0.14 ± 0.15 
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1.64 ± 0.04 ± 0.13 
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0.139 


0.187 ± 0.070 ± 0.097 
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0.048 ± 0.015 ± 0.044 



Table 1 
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